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ABSTRACT 

In analyzing charged particle spectra in space due to galactic cosmic rays (GCR) and solar particle events 
(SPE), the conversion of particle energy spectra into linear energy transfer (LET) distributions is a convenient guide 
in assessing biologically significant components of these spectra. The mapping of LET to energy is triple valued 
and can be defined only on open energy subintervals where the derivative of LET with respect to energy is not zero. 
Presented here is a well-defined numerical procedure which allows for the generation of LET spectra on the open 
energy subintervals that are integrable in spite of their singular nature. The efficiency and accuracy of the numerical 
procedures is demonstrated by providing examples of computed differential and integral LET spectra and their 
equilibrium components for historically large SPEs and 1977 solar minimum GCR environments. Due to the 
biological significance of tissue, all simulations are done with tissue as the target material. 


INTRODUCTION 

In accessing the biological response (risk) to radiation exposure, a customary assumption is that risk is 
related to the energy absorbed per unit mass within a macroscopic volume of the biological material. This implies a 
certain uniformity in the energy deposited to the exposed material, which in general is true for radiation of low 
linear energy transfer (LET) such as x-ray, y-ray, and (3-ray exposures. As experience from exposures with high 
LET radiation was gained, the non-uniformity of the energy deposited in the exposed material was clearly seen as 
an important determinant of biological response, and hence LET dependent relative biological effectiveness (RBE) 
factors came under intense study and provided the basis for the LET dependent quality factor (Q) for use in 
radiation protection (refs. 1-3). As a consequence of such exposure studies, the micro- dosimeter concept was 
developed as an instrument to measure such LET fluctuations. 

Although the concept of an LET dependent RBE is useful for charged particles of relatively low kinetic 
energy, there still are serious concerns about the biological effects of high charge and energy (HZE) ions in deep 
space (outside the earth’s protective geomagnetic field) resulting from the lateral spread of the energy deposited in 
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the exposed material that is larger than a typical cell nucleus (ref. 4). In this case, LET is not sufficient to define an 
RBE, but the track width also plays an important role in calculating the biological response. This was the origin of 
the limited LET concepts (refs. 1 and 5) as related to the lateral spread of the energy deposited in an exposed 
material. The importance of such track structure effects on biological response has been demonstrated for a number 
of biological systems (ref. 6). 

Even though LET as a biological response indicator is most applicable to low charge and low energy 
particles, LET, due to its simplicity, is also considered a rough indicator for HZE particles as well. For this reason, 
we are interested in the generation of differential and integral LET distributions as a guide in identifying 
biologically significant components of HZE. 

In this work, tissue due to its biological significance is used as the only target material. The methodology of 
differential and integral LET simulations as presented here are, however, applicable to any material such as silicon, 
the material of interest in microelectronics. 

The transport of charged ions through tissue is numerically simulated using the deterministic ionizing 
radiation transport code HZETRN (High charge Z and Energy TRaNsport) developed at the Langley Research 
Center (LaRC). This code simulates the radiation shielding response of any material subject to exposure to three 
different environments, namely: the GCR environments in both low earth orbit (LEO) and free space, 
geomagnetically trapped proton fields, and historically large SPEs (ref. 7). 

For GCR at LEO and free space, the environment of solar minimum (1977) is used. This epoch represents the 
worst case scenario (maximum spectral intensities) of the GCR environment. At LEO, the GCR environment is 
further modified (attenuated) by the effects of Earth’s geomagnetic field. The international space station (ISS) 
trajectory (400 kilometer altitude, and 51.6 degrees inclination) is used to simulate the effect of GCR at LEO. The 
ISS case study is chosen here since recent data from the space shuttle flights (space transportation system, STS) to 
the high orbital inclination of ISS show that more than 50% of accumulated dose incurred by ISS crew (except 
during SPE) are from high energy nuclei in the GCR spectra, and of that amount a substantial fraction is from 
highly charged and energetic heavy nuclei. The large contribution of these particles to the accumulated dose is due 
to the high orbital inclination of the ISS at 51.6 degree, since the geomagnetic field strength diminishes near the 
poles, and HZE particles that are deflected at lower latitudes are allowed to reach objects at high inclination LEO 
flights. 

For SPE environments, parametric fits to the February 1956, November 1960, August 1972 (King and LaRC), 
September 29, 1989, and October 1989 events are used (ref. 7). For SPE simulations, only the free space 
environment is considered due to limited ability of SPE particles to penetrate the geomagnetic field except at very 
high latitudes. 


The rest of this report is organized in the following manner: First, a somewhat detailed description of the 
mathematical and numerical modeling of differential LET spectra, its equilibrium component, and stopping power, 
with the connection between the two, as used in this paper is presented. This is followed by a description of the 
integral LET spectra and its equilibrium component. Next, to validate the efficiency of the numerical procedures in 
computing differential and integral LET spectra and their equilibrium components, the results of numerical 
simulations for a number of established SPEs and the 1977 solar minimum GCR environments in free space and 
LEO are presented. Finally, there is a brief concluding remarks section that includes suggestions for future work. 
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DIFFERENTIAL LET SPECTRA 

In radiobiology, relating biological response to the LET of the radiation environment is traditional. For 
example, RBE and the quality factor Q are generally taken as being related to LET. As a consequence, the concept 
of LET spectra has played a role in estimating biological response. Unfortunately, this concept is most useful if the 
flux ( (f) L (L) dL ) of particles with LET (L) between L and L + dL is known. This is generally found by knowing the 
energy flux ( </> e (E) dE ) of particles with energy E between E and E + dE , where L is known as a function of E so 
that and </>(E) are related according to equation 1 below: 


A(L) = 


dL 

dE 


ME) 


( 1 ) 


where the pre-factor 


dL 


dE 


is the Jacobian between the E and L spaces. The difficulty with this approach is that 


dL/dE = 0 at the maxima and minima of the LET curve and that (j) L (L) must be replaced by the sum over the 
various branch functions as: 


A(L) = Z 


dL 

dE 


Me b ) 


( 2 ) 


where E B is the energy of each branch associated with L. That is, for all values of E B the following must hold: 

L = L(E b ) (3) 

Clearly, (j) L (L) does not exist for every value of L but is defined on open intervals not containing values for 
which dL/dE = 0. Furthermore, (j) L (L) is unbounded on the open subintervals over which it is defined, even though 
(j) L (L) is integrable over its domain. From the above arguments, enough challenges obviously exist in finding a 
representation for (p L (L) . This problem can be simplified since L has but one maximum and one minimum other 
than at zero energy. Furthermore, L in the neighborhood of the branch limits has continuous second derivative, 
allowing L to have a Taylor series expansion (approximation) of the form: 

L~L(E h )+ ] -L"(E b )(E-E h ) 2 (4) 

Combining equations 1 and 4, one can show that in the neighborhood of the branch limits </> l (L) 
can be approximated by: 

ML) « ME^lj^a-L^l' (5) 

where the subscript B denotes evaluation at the branch limit. 
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The above considerations are implemented in the following manner. The LET is defined over a numerical grid 
given by the sequence {E. } . The maximum and minimum branch points are found at dL/dE - 0 and are noted by 

E max and ^min > respectively. The sequence {E j } L is defined as those values of E i less than E max , with the main 
sequence {E t } m being defined by E max < E { < E mxn and the sequence {E. } H defined by E mxn < E { . The three 
branch functions are then represented by: 


{*u)b 



E,e{E,} B 


( 6 ) 


where B denotes one of the three branches (that is, B = L, m, or H ). Having a table of values { {(f ) Li , L { ) } B in order to 
reconstruct an adequate representation of the function over each branch may not be sufficient because (j) Li is 
unbounded near the branch limits and an extrapolation into the neighborhood of the branch limit must be provided. 
This is accomplished by recognizing that at the branch point (/> l {L b ) can be approximated by: 


= J \L-L B \~ l (7a) 

and if {E. } is sufficiently close to the branch point, equation 7a can be used to approximate the spectrum by: 

(L ) ~ (/> u ( L b )Jl. — L b | (7b) 

Combining equations 7a and 7b imply that near the branch point (j) L (L) can be approximated by: 

<t>L ( L ) ® 0Li (\ L i ~ L B | /| L - L B |)‘ /2 (70) 

where L is the nearest grid value to the branch limit L B in the appropriate domain. Thus, the data set required to 
reconstruct the LET spectrum is the branch limit values of E m , xx , E mln , L max , and L min and the sequences 
{E t }, {L. }, and {(j) Ei } . Note that the numerical values of the above parameters depend on the charge Z and mass A 
of the particles of the field. Thus, E m , xx , E mln , L max , and L min must be specified for each ion type in the radiation 
field. 


STOPPING POWERS 

As stated earlier, the energy imparted to the medium by a passing charged particle is related to the stopping 
power (the energy loss per unit distance traveled in the medium) and is known as the linear energy transfer (LET), 
where L — S{E ) . To calculate the spectrum of energy deposited, we must evaluate the LET spectrum given by 
equation 1. 

Clearly, a method of generating d(S(E))l d(E) as a function of Eds required. As explained elsewhere (ref. 
8), the stopping powers used in this work are taken from the fitted curves for protons and alpha particles of Ziegler 


4 



at low energy (ref. 9), and are extrapolated using the Bethe formula at high energy. The stopping power of ions with 
a charge greater than 2 is scaled using the effective charge formalism of Barkas (ref. 10), and the corresponding 
alpha particle stopping power of Ziegler. 


In the HZETRN code (ref. 7), the stopping power values are stored in a data array for 15 ion charge values 
ranging from 1 to 92 over an energy grid of 60 points between O.OlMeV and 50 GeV. Furthermore, extrapolation to 
lower energy or higher energy is accomplished by making the stopping power and its first derivative continuous at 


the boundaries of the energy grid. The numerical values are then scaled by 


In 


— yS(lnis,Z) 


to minimize 


numerical errors. The scaled data array is then fitted by a two-dimensional spline with appropriate boundary 
conditions from which d(S(E))l d(E) is found. Figure 1, below shows S(E ) and d(S(E))/ d(E) as functions of 
energy E for a target material of tissue composition for the projectile mass number in the range of 1 to 58, with an 
energy range of 0.01 MeV to 10 GeV . 


S(E) vs. E for GCR Ions in Tissue 



dS(E)/d(E) vs. E for GCR Ions in Tissue 
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Figure 1. Stopping power and its derivative for ion range of 1-58 for a tissue composition. 


The derivatives in the figure above vanish at S(E) maximum at low energy where binding effects dominate 
and S(E) minimum at high energy because of relativistic and density effects. The energies associated with the 
maximum and minimum E max and E min are shown in Figure 2, along with the corresponding values of L max and 

L^n • These provide the basic data for the construction of differential FET spectra. 
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Energy values at branch limits 


LET values at branch limits 



Figure 2. Values of energy and LET at branch points as a function of charge number Z. 


EQUILIBRIUM SPECTRA 

Regardless of whether the source of the projectile ions (GCR, SPE, etc...) is at the boundary or is an interior 
volume source, the accumulated fluence approaches a characteristic equilibrium spectrum that depends on the 
composition of the exposed medium. The equilibrium spectrum for volume sources has been given elsewhere (ref. 
1 1). The fluence within a medium at a distance x from the boundary is given by: 




( 8 ) 


where 0(0, E ) is the fluence at the boundary, E x = R 1 [R(E) + x ] is the energy at depth x, and R(E ) is the range 

of the individual particle with energy E such that E = [/?(£)] . The 0(x, E) term of equation (8) results from 

the formalism (refs. 11, 12) where the balance of high energy gains and low energy losses (ignoring the nuclear 
reactions) is given by: 

a a 


T x ~Te s{e) 


0(x, E) = Nuclear gains and losses 


(9) 


Excluding the nuclear processes, the solution of equation 9 can be expressed as: 

S (E)0(x, E) = S (E x )0(O, E x ) (10a) 

It must be noted that when x » R(E) , then E x ~ /? -1 (x) and a solution to equation (10a) is approximated by: 

(10b) 

S(E ) 

which is the equilibrium spectral limit solution. 

Figure 3, shows the differential LET spectra for the tissue thickness range of 0 < v < 50 cm. for six of the 
historically larger SPEs of solar cycles 19 and 20, where as stated earlier the propagation through the tissue was 
calculated using the HZETRN code. The abundance of high energy protons in the February 1956 spectrum is 
evident with the low LET (high energy) singularity being seen in the figure. Furthermore, note that the existence of 
the higher energy spectra of the LaRC fit to the August 1972 event relative to the King spectrum is evident at the 
largest penetration depths. As is evident in Figure 3, the Bragg peak due to the protons is present at around 1000 
MeV/cm in all six spectrums, while the peak due to the helium secondary production at all depths beyond the 
boundary manifests itself at around 2000 MeV/cm. 
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Freespace dO(L)/d(L) vs. L, for Feb.1956 SPE 



Freespace dO(L)/dL vs. L for Aug. 1972 SPE (King) 



Freespace d<X>(L)/dL vs. L for Sep. 1989 SPE 



Freespace dO(L)/d(L) vs. L for Nov. 1960 SPE 



Freespace dO(L)/dL vs. L for Aug. 1972 SPE (LaRC) 



Freespace dO(L)/d(L) vs. L for Oct. 1989 SPE 



Figure 3. Differential LET spectra for six SPE events of cycle 19 and 20. 
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The differential form of the fit for SPEs as represented by equations 11 a.- Ilf. below, were used as the 
boundary condition to HZETRN to initiate propagation in tissue (ref. 7). 


February 1956: 
November 1960: 
August 1972 (King): 
August 1972 (LaRC): 
September 1989: 


October 1989: 


d(j) P ( £ ) 
d{E) 
d(j) P ( £ ) 
d{E) 
d<(> P ( £ ) 
d(E) 
d(/) P ( E ) 
d{E) 

d(j) p (£) 
d{ E) 


6.0xl0 7 e ' 25 ’ + 9.375x10 V' 320 ' 

£-10 _ ( £— 100 ) 

6.33xl0 8 e 12 + 4.88xl0 6 e 80 

-(— ) 

2.98xl0 8 e 265 


2.20x10' e 


£-100 
7 - 1 30 ) 


+ £(£+ 1876 ) 

1 .446x1 0 8 (—j=======)e 102118 


\/£(£ + 1876) 


(11a.) 


(lib.) 


(He.) 


(lid.) 


£ < 30 MeV (lie.) 


d(j) P (£) 
d(E) 


2.034x1 0 7 
( 938 V 

v t^ + 938 y 


r ^£(£ + 1876) ^ 
v ^30(30 + 1876 ) ; 


£ > 30 MeV 


d(j) p (£) 
d(E) 


6.104x10 8 ( 


£ + 938 
^£(£ + 1876) 


+£(£+ 1876 ) 

92.469 


(Ilf.) 


The corresponding integral form of the fit for SPEs as measured by the ground or space based detectors are 
given by equations 12a-12f. (ref. 7). 


February 1956: 

-(£-1°) - ( £ - 100 ) 

(f) P {> £) = 1.5x1 0 9 e 25 + 3x1 0 8 e 320 


(12a.) 

November 1960: 

£-10 —(£—100) 

(f) P {> £) = 7.6xl0 9 e 12 +3.9xl0 8 e 80 


(12b.) 

August 1972 (King): 

-(— ) 

(f) P {> E) = 7.9xl0 9 e 265 


(12c.) 

August 1972 (LaRC): 

_ ( £- 100 ) 

0 P (>E) = 6.6x1 0 8 e 30 


(12d.) 

September 1989: 

E(E+1816) 

^ P (>£) = 1.477xl0 1 V ' 102118 } 

£ < 30 MeV 

(12e.) 

October 1989: 

(j) P (> £) = 2.94xl0 17 (7£(£ + 1876))" 3 ' 5 

+£(£+1876) 

</> P (> E) = 5.644xl0 10 e 92469 

£ > 30 MeV 

(12f.) 
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Note that the limit of equation 10b, for a pure equilibrium spectrum (any particle type j) at all energies is 
simply given by : 


ME) = 


C 


(13) 


since the numerator in equation 10b becomes a constant for x» R(E) . As an example, the equilibrium 
differential LET spectra for c = 1 is shown in Figure 5 below. Note that the equilibrium spectrum is seen in the high 
LET region for the stopping powers of hydrogen ions (mainly protons) and the helium ions (mainly alpha particles 
produced as nuclear decay products). 


Equilibrium dO(L)/d(L) vs. L spectra for 1977 GCR environment 



Figure 5. Equilibrium differential LET spectra for 1977 GCR solar minimum in tissue. 


INTEGRAL LET SPECTRA 

The integral LET spectrum is given as: 

<£>(> L) = £ Lm “ <f),(L') dL' (14) 

which may be related to integral energy spectra as: 

0(> L) = 0(> E x ) - 0(> E 2 ) + 0(> E 3 ) (15) 

where E v E 2 , and E 3 are the three roots (branch functions) of the relation L = S(E ) . 

The three branch functions of equation (15) are shown in Figure 6 for a range of S(E ) values for a number of 
ions and represent contiguous domains bounded by boundaries where dL! dE = 0. 
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The three branch functions for selected GCR ions in tissue 



Energy (MeV/n) 


Figure 6. The three branch functions for selected GCR ions in tissue. 


The integral LET spectra for the six SPEs discussed previously are shown in Figure 7. To gain perspective, 
we also evaluate the integral equilibrium spectral shape of equation (13) as: 



dE' 

Ee 7 ) 


R^)-R(E) 


(16) 


where we arbitrarily take the variable % to be 10 GeV. The integral LET spectrum for % > E 3 is then: 

0(> L) = R(E 2 ) - R(E X ) + R{£) - R(E 3 ) (17) 

If we examine only the high-LET region, then: 

<f>(>L) = R(E 2 )-R(E l ) (18) 

Results from equation (18) are shown in Figure 8. Thus, the integral LET spectra are characterized by the main 
branch of the LET curve as: 


0(>L)~7?(£ 2 ) 


(19) 


except at the highest LET values where 
C H>L)~(E 2 -E x )IS(E m J 


( 20 ) 


11 




Freespace O(L) vs. L for Feb. 1956 SPE 



Freespace O(L) vs. L for Aug. 1972 SPE (King) 



Freespace <&(L) vs. L for Sep. 1989 SPE 



Freespace O(L) vs. L for for Nov. 1960 SPE 



Freespace O(L) vs. L for Aug. 1972 SPE (LaRC) 



Freespace <X>(L) vs. L for Oct. 1989 SPE 



Figure 7. Integral LET spectra for six SPE events of cycles 19 and 20. 
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Equilibrium ®(L) vs. L spectra for 1977 GCR environment 



Figure 8. Equilibrium integral LET spectra for 1977 GCR solar minimum in tissue. 

As a practical application of the numerical procedures described earlier, the differential LET spectra for the 
1977 solar minimum GCR environment in free space and ISS trajectory are shown in Figure 9 at the tissue depth 
range of 0 < v < 50 cm. As shown, the reduction of high LET components at larger tissue depths increases the 
lower LET component distributions. This is due to the fragmentation of the larger HZE components into smaller 
fragments. Note that there are two types of singularities: the first branch which is discontinuous at higher LET to 
the right of maximum ionization, and the third branch which is discontinuous to the left of minimum ionization. 
The corresponding integral LET spectra for 1977 GCR at free space and ISS trajectory are shown in Figure 10. 


Freespace d<D(L)/d(L) vs. L for 1977 GCR environment iss d<D(L)/d(L) vs. L for 1977 GCR environment 




Figure 9. Differential LET spectra for 1977 GCR solar minimum in tissue. 
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Freespace O(L) vs. L for 1977 GCR environment 


ISS O(L) vs. L for 1977 GCR environment 




Figure 10. Integral LET spectra at 1977 GCR solar minimum in tissue. 


CONCLUDING REMARKS 

The problem of generating linear energy transfer (LET) spectra has been clearly delineated. Although integral 
LET spectra are continuous functions, the differential LET spectra are defined only over open sub-intervals and are 
unbounded near the open end points defined by minimum and maximum LET of specific ion types. An accurate 
method of constructing the spectra over their domain is provided. In future work, the numerical methodologies 
presented here will be used to analyze shielding applications for the protection of biological systems and electronic 
equipments. 
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